Photon spectra from quark generation by WIMPs 



J. A. R. Cembranos, A. Cruz-Dombriz^, A. Dobado, R. Lineros^ and A. L. Maroto 

Departamento de Fisica Teorica I, Universidad Complutense de Madrid, E-28040 Madrid, Spain. 
^ also at Department of Mathematics and Applied Mathematics, University of Cape Town, 7700 South Africa. 
^ IFIC, CSIC-Universitat de Valencia, Ed. Instituts, Apt. 22085, 46071 Valencia, Spain. 

Abstract. If the present dark matter (DM) in the Universe annihilates into Standard Model (SM) particles, it must contribute 
to the gamma ray fluxes that are detected on the Earth. The magnitude of such contribution depends on the particular DM 
candidate, but certain features of these spectra may be analyzed in a model-independent fashion. In this work we provide the 
fitting formula valid for the simulated photon spectra from WIMP annihilation into light quark-anti quark {qq) channels in a 
wide range of WIMP masses. We illustrate our results for the cc channel. 
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I. INTRODUCTION 

According to present observations of large scale struc- 
tures, CMB anisotropies and light nuclei abundances, 
DM cannot be accommodated within the SM of elemen- 
tary particles. Indeed, DM presence is a required com- 
ponent on cosmological scales, but also to provide a sat- 
isfactory description of rotational speeds of galaxies, or- 
bital velocities of galaxies in clusters, gravitational lens- 
ing of background objects by galaxy clusters and the 
temperature distribution of hot gas in galaxies and clus- 
ters of galaxies. The experimental determination of the 
DM nature will require the interplay of collider exper- 
iments and astrophysical observations. These searches 
use to be classified in direct or indirect searches (see [1] 
and references in Introduction in [2]) . Concerning direct 
ones, the elastic scattering of DM particles from nuclei 
should lead directly to observable nuclear recoil signa- 
tures although the weak interactions between DM and 
the standard matter makes DM direct detection extremely 
difficult. 

On the other hand, DM might be detected indirectly, 
by observing their annihilation products into standard 
model particles. Thus, even if WIMPs (Weakly Inter- 
acting Massive Particles) are stable, two of them may 
annihilate into ordinary matter such as quarks, leptons 
and gauge bosons. Their annihilation in different places 
(galactic halo. Sun, etc.) produce cosmic rays to be dis- 
criminated through distinctive signatures from the back- 
ground. After WIMPs annihilation a cascade process oc- 
curs. In the end the stable particles: neutrinos, gamma 
rays, antimatter... may be observed through different de- 
vices. Neutrinos and gamma rays have the advantage of 
maintaining their original direction due to their null elec- 
tric charges. 



This communication precisely focuses on photon pro- 
duction coming from qq channels (except tt channel). 
Photon fluxes in specific DM models are usually ob- 
tained by software packages such as DarkSUSY and mi- 
crOMEGAs based on PYTHIA Monte Carlo event gen- 
erator [3] after having fixed a WIMP mass for the partic- 
ular SUSY model under consideration. In this sense, the 
aim of this investigation is to provide fitting functions for 
the photon spectra corresponding to each individual an- 
nihilation qq channel and, in addition, determine the de- 
pendence of such spectra on the WIMP mass in a model 
independent way. This would allow to apply the results to 
alternative DM candidates for which software packages 
have not been developed. On the other hand, the infor- 
mation about channel contribution and mass dependence 
can be very useful in order to identify gamma-ray signals 
for specific WIMP candidates and may also provide rel- 
evant information about the photon energy distribution 
when qq pairs annihilate. 

Let us remind that the 7-ray flux from the annihilation 
of two WIMPs of mass M into two SM particles coming 
from all possible annihilation channels (labelled by the 
subindex i) is given by: 
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where ((7,v) holds for the thermal averaged annihilation 
cross-section of two WIMPs into two (/''' channel) SM 
particles and p is the DM density. The integral is per- 
formed along the Une of sight (l.o.s.) to the target and 
averaged over the detector soUd angle Ail. 



TABLE 1. bi,ni,ci,d\ and p parameters correspond- 
ing to expression (2) in the cc channel. Mass indepen- 
dent parameters in (2) for this channel are a\ = 5.58 ; 
b2 = 7.90 ; m = 0.686 ; C2 = 0.0 ; 9 = 9.00 • lO"''. 
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III. PROCEDURE 

We have used the particle physics PYTHIA software [3] 
to obtain our results. The WIMP annihilation is usually 
splited into two separated processes: The first describes 
the annihilation of WIMPs and its SM output. The sec- 
ond one considers the evolution of the obtained SM un- 
stable products. Due to the expected velocity dispersion 
of DM, we expect most of the annihilations to happen 
quasi-statically. This fact allows to state that by consid- 
ering different center of mass {CM) energies for the ob- 
tained SM particles pairs from WIMP annihilation pro- 
cess, we are indeed studying different WIMP masses, i.e. 
EcM — 2M. The procedure to obtain the photon spectra 
is thus straightforward: For a given pair of SM particles 
which are produced in the WIMP annihilation, we count 
the number of photons in bins for the variable x = Ey/M. 

Once the PYTHIA simulations have been performed, 
the required parametrization to fit the data for the qq 
channels (except tt) may be written as: 

+ q\n\p{l-x)f^^^^ (2) 

The parameters in expression (2) were considered to be 
WIMP mass dependent. After a fitting process they were 
determined for different WIMP masses, in a range vary- 
ing from 50 to 7000 (or 8000) GeV. Mass dependences 
for the parameters in (2) were fitted by using power laws. 



III. cc CHANNEL 

In order to illustrate the explained procedure, we study 
the cc channel. For this channel there are five mass 
dependent parameters in expression (2): bi,ni, ci, di and 
p presented in Table I. The mass independent parameters 
are ai, ^2. i2> C2 {d2 is thus irrelevant) and q. In Table 



TABLE 2. Fitting power laws in cc channel. 



Parameter 
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Power law(s) 


bi 


50 < M < 8000 


9.90M-*'i3'' 


n\ 


50 < M < 8000 


4.14^-0.148 


c\ 


500 < M < 8000 


0.210M00951 


di 


50 < M < 8000 


1.50M-''-^32 






+ 0.479M-00942 


P 


200 < M < 8000 


8.11M0-812 



II we present the fitting power laws for mass dependent 
parameters. Figure 1 presents spectra for four different 
WIMP masses whereas Figure 2 shows fitting power 
laws for two mass dependent parameters. The results for 
the other qq channels [2] are completely analogous even 
though for each channel, the parameters which are mass 
dependent may be different to the ones for the cc channel. 



IV. CONCLUSIONS 

In this work, we have studied the photon spectra coming 
from WIMP pair annihilation into qq pairs for all the 
charmels (except tt). The covered WIMP mass range was 
from 50 GeV to 8000 GeV. Simulated spectra covered 
the whole accessible energy interval: from extremely 
low energetic photons till photons with one half of the 
available total center of mass energy. 

Once the spectra were simulated, an analytical expres- 
sion (2) was proposed to fit the data. This expression 
depends on both WIMP mass dependent and indepen- 
dent parameters. Our results can both provide a better 
understanding of the DM annihilation channels into pho- 
tons and save an important amount of unnecessary Monte 
Carlo simulations. This fact is particularly important for 
high energy photons, whose production rate is very sup- 
pressed. 

Calculations for all qq channels are available 

at the website http://teorica.fis.ucm. 
es/~PaginaWeb/downloads . html. This work 
was partially supported by MULTIDARK CSD2009- 
00064. 

REFERENCES 

1. K. Sigurdson and M. Kamionkowski, PRL 92, 171302 
(2004); J. A. R. Cembranos et al. , PRL 90, 241301 (2003); 
J. A. R. Cembranos and L. E. Strigari, PRD 77, 123519 
(2008); J. A. R. Cembranos, PRL 102, 141301 (2009). 

2. J. A. R. Cembranos, A. de la Cruz-Dombriz, A. Dobado, 
R. A. Lineros, A. L. Maroto, arXiv: hep-ph/1009.4936. 

3. T. Sjostrand, S. Mrenna and R Skands, JHEP05 (2006) 026 
(LUTP 06-13, FERMILAB-PUB-06-052-CD-T). 



cc channel WIMP annihilation, Myy|,y,p = 100 GeV 



cc channel WIMP annihilation, Muwi^-p = 200 GeV 





CC channel WIMP annihilation, Myyiyp = 1000 GeV 




Total tit - 
Montecarlo simulation 



CC channel WIMP annihilation, Mwimp = 5000 GeV 



0.001 0.01 



Total fit - 
Montecarlo simulation 



FIGURE 1. Photon spectra for four different WIMP masses (100, 200, 1000 and 5000 GeV) in the cc annihilation channel. Red 
dotted points are PYTHIA simulations and solid lines correspond to the proposed fitting functions. 
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FIGURE 2. Mass dependence of c\ and b\ parameters for cc annihilation channel. Crossed points are parameters values found 
after the fitting process for each WIMP mass and solid lines correspond to the proposed fitting functions. 



